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ABSTRACT 

We present the fundamental stellar and planetary properties of the transiting planetary system 
WASP-13 within the framework of the Homogeneous Study of Transiting Systems (HoSTS). HoSTS 
aims to derive the fundamental stellar (T e fi, [Fe/H], M*, i?*), and planetary (Mpi, Rpi, T cq ) physical 
properties of known transiting planets using a consistent methodology and homogeneous high-quality 
dataset. Four spectral analysis techniques are independently applied to a Keck+HIRES spectrum of 
WASP-13 considering two distinct cases: unconstrained parameters, and constrained log g from transit 
light curves. We check the derived stellar temperature against that from a different temperature 
diagnostic based on an INT+IDS H a spectrum. The four unconstrained analyses render results that 
are in good agreement, and provide an improvement of 50% in the precision of T e e, and of 85% in [Fe/H] 
with respect to the WASP-13 discovery paper. The planetary parameters are then derived via the 
Monte-Carlo-Markov-Chain modeling of the radial velocity and light curves, in iteration with stellar 
evolutionary models to derive realistic uncertainties. WASP-13 (1.187±O.O65M ; 1.574±0.048R Q ) 
hosts a Saturn-mass, transiting planet (0.500±0.037M,j up ; 1.407±0.052Rj up ), and is at the end of its 
main-sequence lifetime (4-5.5 Gyr). Our analysis of WASP-13 showcases that both a detailed stellar 
characterization, and transit modeling arc necessary to well determine the fundamental properties 
of planetary systems, which are paramount in identifying and determining empirical relationships 
between transiting planets and their hosts. 
Subject headings: transiting planets — individual: WASP-13 



1. INTRODUCTION 

The detection and characterization of a large number 
of extrasolar planets with a variety of physical proper- 
ties, in different environments and with a range of ages 
is necessary for the understanding of the formation and 
evolution of planetary systems. It is only with precise 
measurements of the fundamental properties of the ex- 
oplanets and their host stars that the planetary bulk 
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composition can be inferred and the planetary structure 
probed, and thus, we can explore the underlying physical 
processes involved in their formation and evolution. 

Tra nsit surveys, such as SupcrWASP (Pollacco et al. 
120061) . have been extremely successful in discovering 
planets for which measurements of their masses and radii 
are possible. These have revealed a large diversity of 
physical properties o f the extrasolar pla nets and their 
host stars (see e.g., iBaraffe et alj |20100 . With more 
than 290 transiting exoplanets confirmed to dat^3, it 
is now possible to conduct statistical studies of plan- 
etary properties, and thus, derive more robust empir- 
ical relationships between the planets and their host 
stars. For example, it is generally thought that in 
the case of Hot Jupiters the planetary radius is cor- 
related with the planet equilibrium temperature and 
the stellar irradiation, and anticorrelated with stellar 



metallicity (e.g . iSantos et al.1 120041: iGuillot et all 120061: 
Laughlin etalT 120111: lEnoch et atl 120121 iFaedi et al.l 



2011UDemorv fc Seaeerl 120111 ). IBuchhave et all (I2012TL 
using recent Kepler results, find that giant planets are 
found around metal-rich stars, while those with smaller 
radii than four times that of the Earth are found to orbit 
stars with a large range in metallicity (-0.6 < [Fe/H] < 
0.5 dex) . This is compatible with previous observational 
results (lUdrv fc Santos! [20071: ISousa et al.ll200l !20Llbt 
IGhezzi et al.l 120 loh . which show that Neptunian planets 
do n ot form preferentia l ly arou nd metal-rich stars. More- 
over, [Adibckvan~^Fj3l] (2012b. a) show that although the 
terrestrial planets can be found in a low-iron regime, they 

14 See http://exoplanct.cu 
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are mostly enhanced by alpha elements as compared to 
stars without detected planets showing that metals con- 
tinue to be important also for the formation of these 
planets. 

This paper presents the pilot study of our project, enti- 
tled Homogeneous Study of Transiting Systems (HoSTS), 
that will derive a homogeneous set of physical proper- 
ties for all transiting planets and their hosts stars, aim- 
ing to minimize the effects of any systematics in the 
measurements due to quality of the data and/or tech- 
nique applied. Individual studies of single planetary 
systems make systematic uncertai nties difficult t o iden - 
tify, and quantify. For example, iMancini et al.l (|201 2T ) 
found HAT-P-8b to have a radius ^14% smaller than 
previous estimates (a difference larger than the quoted 
uncertainties), making it consistent with other transit- 
ing planet radii, and not significantly inflated. As more 
planets are being discovered and/or re-analyzed, the ob- 
served trends, like the anomalously large planetary radii, 
remain to be confirmed. A number of recent studies 
employing consistent analysis procedures for subsets of 
the known st ars with transiting planets have been at- 
temp t ed fe.g.. | Torres et aTll2008t lAmmler-von Eiff et al.l 
[20091: ISouthwortb)l2008t I2009L I2010L I2011L 12012ft : how- 
ever, these have largely employed heterogeneous spec- 
troscopic datasets and adopt the stellar properties, like 
T e g and [Fc/H], from the literature. Because these mea- 
surements are non-homogeneous — arising from different 
spectroscopic analysis techniques applied to spectra ob- 
tained with different spectrographs, different resolution, 
etc. — the typically quoted uncertainties of ~10% in the 
published stellar and planetary mass and radius likely 
contai n currently unchara cterized systematics. More re- 
cently, l^rres^LaD (|2012l ) have thoroughly analyzed new 
and archival spectra (from different instruments, and 
with varied signal-to-noise ratio and resolution) of 56 
transiting planet hosts com paring three differe nt stellar 
characterization methods. iTorres et al.l (|2012l ) focused 
on the stellar hosts, deriving a new set of homogeneous 
spectroscopic stellar properties and have been able to 
identify systematic errors due to the stellar characteriza- 
tion techniques applied. Thus, any empirical trend iden- 
tified among the physical properties, like the observed 
inflated radii of hot Jupiters with respect to planetary 
models, may have to be revised. 

HoSTS extends these previous studies in that the stel- 
lar host properties are derived from a homogeneous, high- 
quality spectral dataset applying four stellar characteri- 
zation techniques, and that the planetary properties are 
also derived consistently By means of our homogeneous 
spectral dataset and subsequent analyses, we will be able 
to investigate systematic uncertainties on the derived 
stellar properties aris ing not only from the methodol- 
ogy, as exemplified bv lTorres et al.l ()2012[ ). but also from 
the quality of the data. We will combine iteratively our 
results with the best available radial velocity data and 
transit photometry in the literature to derive a homoge- 
neous set of properties for the transiting systems. The 
resulting consistent set of physical properties will allow 
us to further explore known correlations, e.g., core-size of 
the planet and stellar mctallicity, and to newly identify 
subtle relationships providing insight into our fundamen- 
tal understanding of planetary formation, structure, and 
evolution. And thus, this will allow us to reevaluate the 



planetary properties, of each planet alone, and with re- 
spect to different planet populations. 

In this paper, we present our HoSTS pilot study of 
the transiting system WASP-13 which is composed of a 
Saturn-mass planet around an early G-type star. Section 
[5] describes both the data acquired by our team and the 
data from the literature utilized in our analyses. Sec- 
tion [3] describes the spectral analysis of the WASP-13 
spectra implementing four different techniques in order 
to derive the stellar spectroscopic properties. For each 
of the four stellar characterization methods, we present 
two different cases for which stellar properties have been 
derived: a) an unconstrained analysis where all parame- 
ters are left free; and b) applying an external constraint 
on log g. We have used the temperature diagnostic based 
on H a to verify the values derived from the stellar char- 
acterization methods, as well as the effect of fixing T c g 
on the other spectroscopically determined stellar prop- 
erties. Then we describe the modeling of the system's 
radial velocity and light curves to derive the stellar and 
planetary properties. In section [4] we discuss the results 
from our pilot study of WASP-13, and outline the future 
work of the HoSTS project. 

2. DATA 

For our analysis of the planetary system WASP-13, 
we have acquired new spectroscopic data described be- 
low in £12.21 and £12.11 and have utilized the data avail- 
able from the literature and the SupcrWASP archive. 
The SupcrWASP light curve includes over 12,100 data 
points observed with the SuperWAS P-North facility in 
La Palma, Spain (jPollacco et al.ll2006l ) from 2006 Novem- 
ber to 2009 April (see Fig. [TJ left panel), These data 
expand the span of the SuperWASP light curve from 
that of the discovery paper for two additional years, and 
have a median pho t omet ric uncertainty of 0.006 mag. 
From iSkillen et all (pOOl ). we have obtained the sys- 
tem's radial velocities acquired with the SOPHIE in- 
strument mounted on the 1.9-m telescope at the Haute 
Provence Observatory (see Fig. UJ right panel), as well as 
the James Gregory Telescope (JGT) differential photom- 
etry in the i?-band. Additionally, we have adopted the 
high-cadence, high-precision light curves observed with 
the RIS E instrument on the Liverpool Telescope pub- 
lished bv lBarros etaLl (|2012t ). All follow-up light curves 
are shown in Fig. [2] 

2.1. WASP-13 H a Long-slit Spectrum 

For this work, we obtained a long-slit spectrum of the 
planet host WASP-13 around the H a line using the In- 
termediate Dispersion Spectrograph (IDS) mounted on 
the 2.5-m Isaac Newton Telescope at the Roque de los 
Muchachos Observatory in La Palma, Spain. We used 
the H1800V grating with the RED+2 CCD, and a T/4 
slit yielding a dispersion of 0.35 A/pixel and a resolution 
of R-10 000 at 6560A. The observations of WASP-13 and 
the standard calibrations, including a spectrum of Vega 
and Arcturus, were taken on 2012 May 08 and 09. To 
allow for a precise measurement of the H a profile, the 
WASP-13 spectrum had a signal-to-noise ratio of ~500 
as calculated by the IDL function DER_SNR. The derived 
temperature f ^3.1[) is used as an independent check on 
the effective temperature of WASP-13 derived from the 
stellar characterization methods. 
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Fig. 1. — Left: SuperWASP-N Light Curve of Transiting Planetary System WASP-13. Comprising over 12,100 data points and s pann ing 
from 2006 November to 2009 April, the SuperWASP light curve is shown with gray points with the model transit light curve (see £13-31) in 
black. The photometric data has been folded over the orbital period, and exhibits at phase zero the characteristic dip in brightness as the 
Saturn-mass planet WASP-13b [M„\ = 0.50 ± 0.01 Mj up ) passes in front of its host star blocking part of the stellar light every 4.353 days. 
The median photometric uncertainty of the light curve is .006 mag. Right: The radial velocity measurements from the SOPHIE instrument 
at the Observatoire d'Hautc Provence (Skillen ct al. 2009) arc shown against the our model radial velocity curve that describes the reflex 
motion of WASP-13 due to the presence of the planet. The residuals to the fit are shown in the lower panel. 
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Fig. 2.— WASP-13b Follow-up Transit Light Curves. The transi t 
light curve s of the WASP-13 s ystem from JGT HSkillen et al.ll2T)09T) 
and RISE fBarros ct al. 20l3) are shown in red points overplotted 
with the mode l light curves (in solid black lines) from our MCMC 
analysis f ij3.3t . with their corresponding residuals including the 
uncertainty of the photometric data directly below (red errorbars). 
The light curves and residuals have been shifted vertically from 
zero for clarity. 



2.2. WASP-13 HIRES Spectrum 

We observed WASP-13 on 2011 March 14 UT with the 
H igh Resolution Ec helle Spectrograph (HIRES) on Keck- 
I (jVogt et al.lll994D . We observed in the spectrograph's 
"red" (HIRESr) configuration with an echelle angle of 
-0?018 and a cross-disperscr angle of 0?737. We used 
the KV418 order-blocking filter and the 0'.'57x7'.'0 slit, 



and the chip was binned by 2 pixels in the spatial di- 
rection during readout. The resulting resolving power is 
R-72,000. 

We obtained three consecutive integrations of WASP- 
13, each of 600 s. ThAr arc lamp calibration exposures 
were obtained before and after the WASP-13 exposures, 
and sequences of bias and dome flat-field exposures were 
obtained at the end of the night. The WASP-13 expo- 
sures were processed along with these calibrations us- 
ing standard IRAP0 tasks and the MAKEE reduction 
package written for HIRES by T. Barlow. The latter 
includes optimal extraction of the orders as well as sub- 
traction of the adjacent sky background. The three ex- 
posures of WASP-13 were processed separately and then 
median combined with cosmic-ray rejection into a sin- 
gle final spectrum. The signal-to-noise ratio of the final 
spectrum is ^300 per resolution element. 

3. ANALYSIS 

In this section, we describe the methods applied to our 
WASP-13 dataset in order to derive the physical proper- 
ties of the planetary system. 

3.1. Determination o/T e g from H a Spectrum 

This analysis is based on the H a spectrum of WASP- 
13 described in ^2.11 The Balmer lines provide an excel- 
lent T e g diagnostic for stars cooler than about 8 000 K 
due to their virtually nil gravity dependence dGravl 
l2008h . Normalization of the observations is critical, 
which the shape of the B almer line must be preserved 
(jSmith fc Dwor ctskv 198ct). The extracted spectrum was 
normalized using a low-order polynomial fitted to the 
continuum regions more than 100 A either side of H a , 
in order to avoid any distortion due to the weak wings 
of this profile. The spectrum was then analysed using 

15 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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uclsyn with lCastelli et al.l (|1997l) ATLAS9 models with 
no over shooting and H a profiles calculated using VCS 
theory (|Vidal et al.l 119731 ). The best- fitting H a profile 
has T eff = 5850±60 K. 

However, the use of Balmer lines as temperature 
diagnostics is not without its difficulties due to un- 
certain ties caused by different line broadenin g the- 
ories dStehle fc Hutcheonl H99l iBarklem et all 120001; 
lAllard et al. 120081). the treatment of atmo s pheri c con- 
vection (|Gardiner et al] 119991: iHeiter et al.l l200l . and 
non-LTE effects (jBarkleml |2007D ~ We, therefore, also 
fitted the Hn pro file in the KPNO solar spectrum 
(jKurucz et al.fll984l) which gave a T oS 70 ± 20 K lower 
than the direct value of 5777 K. Thus, H a appears to 
underestimate stellar effective temperatures. Adding 70 
± 20 K to the T of f derived above, it gives T cS = 5920 ± 
60 K f or WASP- 13. This w as recently investigated in de- 
tail bv lCavrel et al.1 (|2011| ) who found similar systematic 
differences, and provide a correction: 

T c ff (direct) = 20.3 + 1.014 x T cS (H a ) 

with an uncertainty of 31 K. Applying this correction, 
and adding the errors in quadrature, gives T c ft = 5950 ± 
70 K for WASP-13. 

Both temperatures from our H a temperature diag- 
nostic agree with each other, and with the T e g derived 
from the infrared flux method (T c g = 5935 ± 183 K; 
iSkillen et al.l l2009h . Moreover, they are consistent with 
the temperatures derived in the stellar characterization 
methods described below. 

3.2. Stellar Characterization Analysis 

We apply to the HIRES echelle spectrum ( $12.2$ four 
different stellar characterization methods that are exten- 
sively used in the exoplanet literature. Each method is 
done independently from each other and is described in 
the subsections below ( 33.2.1ff3T2l| . Method A is based 
on the technique of spectral synthesis, which compares 
an observed spectrum to synthetic model spectra gener- 
ated for a range of stellar parameters. The best fitting 
model (based on a \ 2 minimization) defines the final at- 
mospheric parameters. The other three methods (B, C, 
and D) are based on the principle of excitation/ionization 
equilibrium of iron lines, in which equivalent width mea- 
surements of many lines are used to determine iron abun- 
dance, and the stellar atmospheric properties. Meth- 
ods B, C, and D are each unique in their choice of 
linclists, model atmospheres, equivalent width measure- 
ments, continuum normalization, and convergence cri- 
teria. Furthermore, Method B makes an absolute iron 
abundance measurement of the star, whereas the other 
three methods are differential analyses, and derive a stel- 
lar metallicity relative to the Sun. Method A includes a 
careful determination of the linelist parameters, i.e., ex- 
citation potential and oscillator strength, to match the 
spectrum of the Sun using spectral synthesis. Method 
C docs a differential line-by-line analysis relative to the 
Sun using the same instrument setup. While, Method 
D uses measured equivalent widths and a standard Solar 
iron abundance (e.g., 7.50 ± 0.04: lAsplund et al.ll2009l) 
to derive the line properties. 

For each of the four stellar characterization methods, 
we present two distinct cases: a) the unconstrained anal- 



ysis, where T e ff, logg, and [Fc/H] (and vshii when ap- 
propriate) arc derived freely; and b) constraining logg 
from the mean stellar density as determined from the 
transit model to derive the other stellar properties. The 
spectroscopically-determined stellar parameters derived 
for each case above with all four methods are shown in 
Table [1] Furthermore, the T e g derived using a differ- 
ent temperature diagnostic based on an H a spectrum is 
used to check the stellar properties in Table Q] as well 
as to explore the effect on the derived stellar spectro- 
scopic properties by fixing T e g. The resulting spectro- 
scopic properties are shown in Table [5] 

Given the quality of our data, and the nuances of 
each of the stellar characterization methods, the pre- 
ferred solution for each of the four methods is that de- 
rived through the unconstrained analysis. In the last 
column of Table [TJ we present the weighted mean of 
the four unconstrained solutions for each method which 
are used to derive the stellar mass and radius, and the 
planetary properties (see £|3.3[) . Additionally, we report 
two errors on the stellar properties of Table [1] the first 
is calculated from the weighted quadrature sum of the 
individual internal errors, and the second is a measure 
of the systematic uncertainty based on standard devi- 
ation of the individual measurements. The systematic 
uncertainty is likely to be underestimated in the case of 
WASP-13, because only the four measurements from the 
unconstrained cases are taken into account. A more real- 
istic systematic uncertainty will be possible once a larger 
HoSTS sample has been analyzed in the same consistent 
manner as we present in this paper. 

3.2.1. Method A - SME 

Method A consists on the imple mentation of Spec- 
trosc opy Made Easy (SME v. 3. 54 iValenti fc Piskunovl 
I1996T ) to derive stellar parameters of WASP-13 described 
below. We base the g eneral method of o ur SM E anal- 
ysis on that given in IValenti fc Fischerl (120051 ) includ- 
ing the grid of model atmospheres and derivation of 
macroturbulcnce; however, we use a line list, synthesized 
wavelength ran g es, an d abundance patter n adapted from 
IStempels et all ((20071) : iHebFeFall (120091) . 

In general, SME uses the Levenberg-Marquardt (LM) 
algorithm to solve the nonlinear least-squares problem 
of fitting an observed spectrum with a synthetic spec- 
trum. Like any nonlinear least-squares algorithm, the 
LM based solver in SME requires a good initial guess 
and a smoothly varying \ 2 surface in order to consis- 
tently find the absolute global minimum (what we are 
calling the optimal solution). In addition, a single SME 
best-fit solution does not allow for an estimation of the 
error in the solution apart from the error calculated from 
that solution's covariance matrix. This does not take into 
account the deviations from the best-fit solution depend- 
ing on the specific choice of initial parameter values nor 
the internal precision of the solver. 

We have expanded on the technique outlined in 
IValenti fc Fischerl (|2005D that allows us to operate SME 
in an automated fashion and explore the effect of the 
initial conditions on the final resulting stellar parame- 
ters. Using the ACCRE High-Pcrformancc Computing 
Center at Vandcrbilt University, we have developed an 
extensive Monte Carlo approach to using SME. We start 
by randomly selecting 500 initial parameter values from a 
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multivariate normal distribution with 5 parameters: T e g, 
logs, [Fe/H], [M/H], and vsini. For WASP-13, we de- 
fined this distribution usin g the derived stella r parame- 
ters and uncertainties from lSkillen et all ()2009f ). The mi- 
croturbulence (v t ) for each of these initial values in this 
multivariate distribution was fixed at 1.01 ± 0.17 kms -1 , 
estimated using a polynomial fit (Eq. QJ to the HARPS 
sample of stellar T e g and microturbule nce (jSousa et al.l 
1201 lbf) at T eff from lSkillenetaLl (|2009t) . The value of v t 
is kept fixed throughout our SME analysis, because the 
scatter of the HARPS sample around this temperature 
is larger than the change in vt in the range of tempera- 
tures explored in all cases. Furthermore, the change in 
the microturbulence value within its uncertainties does 
not affect significantly the derived [Fe/H]. 

v t = 0.909148 + (T cff -5700)/1318 + (T cff -5700) 2 /1660 2 

(1) 

We then allow SME to find a best-fit synthetic spec- 
trum and solve for the free parameters for the full distri- 
bution of initial guesses, producing 500 best-fit solutions 
for the stellar parameters. We determine our final mea- 
sured stellar properties by identifying the output param- 
eters that give the optimal SME solution (i.e., the solu- 
tion with the lowest x 2 ). The overall SME measurement 
uncertainties in the final parameters are calculated by 
adding in quadrature: 1) the internal error determined 
from the 68.3% confidence region in the \ 2 map, and 2) 
the median absolute deviation of the parameters from the 
500 output SME solutions to account for the correlation 
between the initial guess and the final fit. 

Following this procedure, we solved for the parameters 
of WASP-13 first letting all fitted parameters be free; 
then using a constraint on log g, and letting the other pa- 
rameters free; and, finally analyzing the spectrum fixing 
T c ff . The resulting optimal parameters and uncertainties 
are given in the first results column of Tables [1] and [3J 

Our choice of the unconstrained analysis as the pre- 
ferred solu tion for method A does not follow the con- 
clusions of iTorres et al.l (|2012f) . where they suggest that 
fixing the log g (case b above) is the best approach when 
using synthesis-based methods. However, we find that 
our method A differs in their implementation of SME in 
the treatment of microturbulence, the linelist, the sam- 
pling of a large parameter space in initial parameters, 
and the convergence criteria. Specifically our linelist in- 
cludes the Na I D region between 5849-5950 A, a gravity 
and temperature sensitive line, as well as the gravity sen- 
sit ive Mg b triplet reg ion. In similarity to the ana lyses 
by ITorres etafl (|20iaf and lValenti fe Fischer! (|200l . we 
did not include the H a region due to the difficulty in 
normalizing the continuum for such a broad line in an 
echelle spectrum. A more in-depth comparison will be 
possible on the larger HoSTS sample, with which we will 
be able to state more robustly whether we need to con- 
strain log g or not, as well as to explore the dependence 
of physical properties on the data. 

3.2.2. Method B - UCLSYN 

Method B consists in the analysis performed with 
the spectral synthesis package UCLSYN (Universi t y Col- 
lege L o ndon SYNthesis:lSmith fc Dworetsky|1988tlSmithl 
119921: iSmallev. Smith, fc Dworetskvl 120011 ) using the 



methods given in lDovle et al.l (|2013l ). In general, the sur- 
face gravity (logg) was determined using the ionisation 
balance of the Fe I and Fe II lines, as well as from the Ca I 
line at 6439A and the Na I D lines. The excitation bal- 
ance of the Fe I lines was used to determine the effective 
temperature (T B s). A null dependence was required be- 
tween Fe abundance and equivalent width in order to as - 
certain the microturbulence (vt) using the lMagainl(|1984l ) 
method. The Fe abundance was determined from equiva- 
lent width measurements of several unblended lines, and 
additional least squares fitting of lines was performed 
when required. The projected stellar rotation velocity 
(v sin i) was determined by fitting the profiles of several 
unblended Fe I lines. A value for macroturbulence (v mac ) 
of 3.0 ±0.3 kms" 1 was as sumed, based on the calibra- 
tion bv lBruntt et al.l ([20101 ). 

The parameters obtained from the analysis are listed 
in the second column of Tables Q] and [2] 

3.2.3. Method C - ARES/MOOG + Schuler linelist 

In the case of the Method C for the unconstrained anal- 
ysis, the atmospheric parameters (T e ff, log g and Vt) and 
metallicity ([Fe/H0) of WASP-13 were derived using the 
standard spectroscopic method based on the excitation 
and ionization equilibrium of Fe I and Fe II lines. The 
[Fe/H] abundances were normaliz ed on a line-by - line ba - 
sis to the solar values taken from lSchuler et aLl (|2011al) . 
The analysis was done in Local Thermodynamic Equilib- 
rium (LTE) using the 2010 version of MOOGO (|Snedenl 
[197! and one-dimensional plane-parallel model atmo- 
spheres interpolated from the O VER grid of ATLAS9 
models (|CasteIli fc KuruczH2004f ). 

The line list was adopted from ISchuler et all (I2011bfl 
and the equivalent widths ( EWs) were measur ed using 
the automatic code ARES (jSousa et al J 120071) . Effec- 
tive temperatures and microturbulence velocities were 
iterated until the slopes of [Fc/H] versus \ (the excita- 
tion potential of the lines) and log(£W/A) (their reduced 
equivalent widths) were respectively zero; i.e., until the 
individual [Fc/H] abundances were independent of ex- 
citation potential and reduced EWs. Surface gravities 
were iterated until the [Fe/H] abundances determined 
from Fe I and Fc II lines were equal. The iteration of the 
atmospheric param eters was done auto matically, using 
codes adapted from lGhezzi et ail (|2010l ). Any lines with 
[Fe/H] abundances that deviated more than 2a from the 
average were removed and the above iteration was re- 
peated until convergence was achieved. The final line 
list contained 45 Fe I and 5 Fe II lines. 

The internal uncertainties on the atmospheric parame- 
ters were estimated as follows. The error of the microtur- 
bulence was determined by varying this parameter until 
the slope of [Fe/H] versus \og(EW/\) was equal to its 
standard deviation. The uncertainty of the effective tem- 
perature was obtained by changing this parameter until 
the slope of [Fc/H] versus x was equal to its standard 
deviation. The error of vt was also taken into account 
when calculating the uncertainty of T c g . The error of the 
surface gravity was obtained by varying this parameter 
until the difference between the average [Fe I/H] and [Fc 

16 [Fe/H] = A(Fe I)* - A(Fe l) , where A(Fe I) = log 
[N(Fe I)/N(H)] + 12 

17 Available at http://www.as.utcxas.edu/~chris/moog.html 
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TABLE l 

Spectroscopically-determined Stellar Parameters of WASP-13 



Weighted Mean 



Unconstrained 



T cB (K) 
logs 

log A(Fc) 
[Fe/H] 

usini (kms 1 ) 
vt (kms -1 ) 



6003 ± 65 
4.16 ± 0.08 
7.54 ± 0.06* 
0.04 ± 0.05 
5.79 ± 0.08 
1.01 ± 0.17* 



5955 ± 75 
4.13 ± 0.11 
7.60 ± 0.09 
0.10 ± 0.09+ 
5.26 ± 0.25 
0.95 ± 0.10 



5919 ± 30 
4.02 ± 0.06 
7.54 ± 0.05t 
0.04 ± 0.02 



6025 ± 21 
4.19 ± 0.03 
7.58 ± 0.05t 
0.08 ± 0.02 



1.53 ± 0.09 1.28 ± 0.10 



5989 ± 16 ± 48 
4.16 ± 0.03 ± 0.07 
7.56 ± 0.03 ± 0.03 
0.06 ± 0.01 ± 0.03 
5.74 ± 0.08 ± 0.38 
1.27 ± 0.06 ± 0.29 



Fixing logs = 4.10 ± 0.04 dex 



T ofi (K) 
log A(Fe) 
[Fe/H] 
tjsini (kms 
vt (kms -1 ) 



5994 ± 150 
7.55 ± 0.12t 
0.05 ± 0.11 
5.86 ± 0.22 
1.01 ± 0.17* 



5955 ± 70 
7.59 ± 0.09 
0.09 ± 0.09* 
5.26 ± 0.25 
1.00 ± 0.10 



5912 ± 30 
7.54 ± 0.05* 
0.04 ± 0.02 



6048 ± 63 
7.56 ± 0.07* 
0.06 ± 0.06 



1.53 ± 0.09 1.33 ± 0.10 



* The first error {a w ) is derived from the uncertainties in the individual measurements (o"j), 
&w = [ 53 l/o"i J -1 / 2 , and the second error is calculated from the standard deviation of the 
individual measurements. 

* Deri ved using the most c urrent value for the solar abundance of iron, Alog(Fe)© = 7.50 ± 
0.04 HAsplund et al.|[2009l1 . 

* The vt has been adopted from the empirical relationship described in £13.2,11 and is not 
included in the weighted mean. 



TABLE 2 

Stellar Properties of WASP-13: Fixing T eff 



5950 ± 70 K 



B 



C 



1) 



logs 

logA(Fe) 
[Fe/H] 

v sin i (km s ) 
vt (kms -1 ) 



4.14 ± 0.19 
7.56 ± 0.07* 
0.06 ± 0.06 
5.88 ± 0.05 
1.01 ± 0.17 



4.13 ± 0.11 
7.60 ± 0.09 
0.10 ± 0.09* 
5.26 ± 0.25 
0.95 ± 0.10 



4.07 ± 0.13 
7.56 ± 0.06* 
0.06 ± 0.04 



4.06 ± 0.10 
7.52 ± 0.07* 
0.02 ± 0.06 



1.53 ± 0.08 1.30 ± 0.10 



* Derived using the most c urrent value for the s olar abundance of iron, 
Alog(Fe) Q = 7.50 ± 0.04 iTAsplund ct al.11200911 . 



II/H] abundances were equal to the standard deviation 
of the latter (divided by the square root of the number of 
Fe II lines) . The contribution from T e g was also included. 
Finally, the uncertainty of [Fe/H] is a combination of the 
standard deviation of the [Fc I/H] abundance (divided 
by the square root of the number of Fe I lines) and the 
variations caused by the errors in T e ff, logg, and Vt, all 
added in quadrature. We note that these are the internal 
errors of the spectroscopic differential analysis used here 
and that the real uncertainties (e.g., from the comparison 
with other similar studies) might be larger. 

In the second case, the surface gravity is fixed to the 
value determined from the analysis of the mean stellar 
density (logg = 4.10 ± 0.04), the other parameters are 
iterated upon (with the same line list as above). The 
uncertainties of T ff and u t were estimated as described 
above. The error of the metallicity took into account 
all four contributions described above, but the influence 
of the surface gravity was estimated by varying log g by 
± la (i.e., fixing this parameter at the values 4.06 and 



4.14) and iterating the atmospheric parameters again. 
The larger difference between the new value and the one 
obtained with logg = 4.10 was taken as the error on 
[Fe/H] due to the uncertainty on the surface gravity. We 
have observed that this variation of log g has no effect on 
[Fe/H] . 

In the third case, we fixed the effective temperature 
to the value determined from the analysis of the H a line 
(T c ff = 5950 ± 70 K), and iterated the other parameters 
(with the same line list as above). The uncertainty of 
Vt was estimated as in the original procedure (with free 
parameters) . The influence of the effective temperature 
on the errors of the surface gravity and metallicity was 
determined by varying T e g by ± la (i.e., fixing this pa- 
rameter at the values 5880 and 6020 K) and iterating the 
atmospheric parameters again. The larger differences be- 
tween the new values and the ones obtained with T e g = 
5950 K were taken as the errors on [Fe/H] and logg due 
to the uncertainty on the effective temperature and then 
added in quadrature to the other contributions, giving 



HoSTS : WASP- 13 



7 



the final values quoted above. 

3.2.4. Method D - MOOG/ARES + Sousa linelist 

In the case of Method D, the spectroscopic parame- 
ters were derived starting with the automatic measure- 
ment o f equivalent width s of Fe-I and Fe-II lines with 
ARES {S 

ousa et al.1 12007ft and then imposing excitation 
and ionization equilibrium using a spectro scopic analysis 
in LTE with the help of the code MOOG (jSnedenlll973ft 
and a GRI D of Kurucz A tlas 9 plane-parallel model at- 
mospheres (|Kuruc 

The Fe I and Fe II line list is composed of more than 
300 lines that were individually tested in high resolution 
spectra to check it s stability to an a utomatic measure- 
ment with ARES (|Sousa et all 12008ft . The atomic data 
of the line s were obtained fro m the Vienna Atomic Line 
Database (jKupka et al.lll99"9ft but the oscillator strength 
(log gf) of the lines were recomputed through an inverse 
analysis of the Solar spectrum allowing this way to per- 
form a differential analysis relatively to th e Sun. A full 
descri ption of the method can be found in iSantos et al.l 
(|200l and lSousa et~aTl (12008ft . 

We have only reported in Tables Q] and [2] the internal 
errors derived from our method. Typically, we report 
a more realisitic uncertainty that considers the typical 
dispersion pl otted in each comp arison of parameters, as 
presented in iSousa et alj (|2008ft . A more complete dis- 
cussion about the systematic errors generall y derived for 
this sp ectroscopic method can be found in ISousa et al.l 
(|2011aft . However, given that in this paper we derive the 
systematic uncertainty from the comparison against the 
resulting spcctroscopically-determined parameters from 
Methods A, B and C, we only report the internal errors 
for each of the three cases. 

For the constrained cases, the same method was used 
but fixing each specific parameter (logg and T e ff) in the 
process. We also used the same procedure for the deter- 
mination of the errors for which the uncertainty in each 
constrained parameter was considered. 

3.3. Transit Model 

The planetary properties were determined using a si- 
multaneous Markov-Chain Monte Carlo (MCMC) analy- 
sis including the WASP photometry, and the high preci- 
sion photometry, together with the radial velocity mea- 
sureme nts. A detailed descr i ption of the method is 

;iven in lCollier Cameron et al.l (|2007ft and lPollacco et al] 

2008ft . 

Our iterative fitting method uses the following param- 
eters: the epoch of mid-transit To, the orbital period P, 
the fractional change of flux proportional to the ratio of 
stellar to planet surface areas AF = R^/R^, the transit 
duration T14, the impact parameter b, the radial velocity 
semi-amplitude K\ , the stellar effective temperature T c g 
and metallicity [Fe/H], the Lagrangian elements 1/ecos w 
and ^/esinw (where e is the eccentricity and u> the lon- 
gitude of periastron), and the systematic offset velocity 
7. The spectroscopically determined T c s and [Fc/H] pre- 
sented in the last column of TablcQ]with the two reported 
errors added in quadratures are used within our MCMC 
code as priors in the stellar mass determination from the 
empirical Torres/Enoch relationship (see below). The 
sum of the x 2 f° r au input data curves with respect to 
the models was used as the goodness-of-fit statistic. 



An initial MCMC solution with eccentricity as a 
free parameter, was explored for WASP-13 deriving a 
small eccentricity (e = 0.10). However, the probabil- 
ity that it is a spurious non-circular orbit as defined by 
ILucv fe Sweenevi (|1971l ) is 1.0 in agreement with previ- 
ous circular solutio ns for the system (jSkillen et al.ll2009l : 
iBarros et al.ll2012"ft . Thus, we adopt a circular orbit for 
the rest of our analysis. For the treatment of the stellar 
limb-darkening, the four-coefficient law of iClaretl (|2000t 
12001 was used with their derived coefficients in the R- 
band for both the JGT and as an approximation for the 
WASP photometry, which is in a V + R passband. In 
the case of the RISE photometry, we used the same limb- 
darkening law, using the coefficients derived specifically 
for the RISE passband and CCD response by I. Howarth, 
following the procedure in lHowarthl ()2011h . 

From the parameters mentioned above, we calculate 
the mass M, radius R, density p, and surface gravity 
log g of the star (which we denote with subscript *) and 
the planet (which we denote with subscript p i), as well 
as the equilibrium temperature of the planet assuming 
it to be a black-body (T p i i a=o) and that energy is ef- 
ficiently redistributed from the planet's day-side to its 
night-side. We also calculate the transit ingress (egress) 
times 712(134), and the orbital semi-major axis a. These 
calculated values and their 1-a uncertainties from our 
MCMC analysis are presented in Table [3] The observed 
light curves are plotted against the model light curves 
with their residuals in Fig. [2] 

The stellar mass of planet host star has been de- 
rived within our M CMC analysis from the empirical 
iTorres et al.l (|2010ft calibration, which is based on the 
precisely measured masses and radii of eclipsing binary 
stars, and relates log<?, [Fe/H], and T c g to the stellar 
mass and radius. However, while T e g can be determined 
with high precision from the stellar spectrum (see §3.11 
and £|3.2p . log g is usually poorly constrained, and thus 
stellar masses derived from the spectroscopic log g can 
have large uncertainties and can suffer from systematics. 
Thus, our MCMC method derives the stel lar mass us- 
ing th e empirical calibration as descr ibed bvlEnoch et al.1 
(|2010ft . which is the based on that bv lTorres et al.1 (12010ft 
but relies on the directly measured p+, instead of logg. 
The stellar density, p*, is directly determined from tran- 
sit light curves and as such is independent of the stellar 
mass, and the effective temperature determined from the 
spectrum (|Sozzetti ct al. 200"7t lHebb ct al. 200j|), as well 
as of theoretical stellar models (M p i <C M* is assumed; 
see iSeaeer fe Mallen-Ornelasl 12003ft . The error on the 
stellar mass that we first derived from the MCMC, which 
is based on the empirical relationship, seemed underes- 
timated (~1%), and thus so are the errors of the other 
properties that depend on the stellar mass (e.g., orbital 
separation). Therefore, we adopt a more realistic uncer- 
tainty in the stellar mass from the comparison of stellar 
evolutionary models to the observed properties of WASP- 
13, as described in the paragraph below and is shown in 
Fig.reftracks. This uncertainty in the stellar mass is in- 
cluded in the final MCMC analysis (see Table [3]) and is 
propagated through all other dependent parameters. 

Because the planet physical properties depend directly 
on the stellar ones and to assure the validity of our 
MCMC results, we have derived independently the stel- 
lar mass from the Yonsei-Yalc (Y 2 ) stellar evolutionary 
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models (see Fig. 131 iDemarque et al"]|2004D . Furthermore, 
this consistency check allows us to estimate the age of 
the planetary system, and to derive a realistic error on 
the mass of the stellar host. We have interpolated the Y 2 
models considering the 1— a errors in the measured p*, 
and in the spectroscopically determined [Fc/H] and T e g. 
As shown in Fig. [31 the evolutionary state of WASP- 
13 is not uniquely determined by the measured stellar 
properties alone. Different mass tracks for the pre-main 
sequence, main sequence and post-main sequence evo- 
lutionary phases overlap in the stellar density-effective 
temperature-metallicity plane at the position of WASP- 
13. We must use additional criteria to identify the most 
likely mass and age for WASP-13. First, there is no 
evidence of youth in our WASP-13 data; the measured 
lithium abundance, A(Li) = 2. 11 ± 0.08 dex, is consisten t 
with an age of several Gyr (see lSestito fc Randich 20051 ). 
Thus, we do not consider the pre-main sequence phase as 
a plausible evolutionary state for WASP-13. In addition, 
the measured surface gravity and temperature rule out 
the post-main sequence, red giant phase. Therefore, the 
most likely scenario is that WASP-13 is at the end of its 
main sequence lifetime and may or may not have reached 
the phase of overall contraction before exhausting hydro- 
gen in its core. According to these models, the stellar 
mass is between 1.175, and 1.245 M Q depending on the 
precise phase of evolution, and the age of the system is 
between 4-5.5 Gyr. Thus, the stellar mass derived from 
the Y 2 theoretical models is consistent with the stellar 
mass from our MCMC analysis (1.187 Mq). However, 
the range of possible stellar masses derived from a single 
set of evolutionary models is larger than the uncertainty 
on the stellar mass given by the empirical Enoch relation. 
Therefore, we conservatively adopt a larger uncertainty 
of cr^ = ±0.065 M© on this parameter to account for 
all plausible mass values. 

4. DISCUSSION AND CONCLUSIONS 

In our re-analysis of WASP-13 presented in this paper, 
we are able to more accurately determine the system's 
physical properties than in previous studies. An improve- 
ment of 85% and 50% in stellar metallicity and effec- 
tive temperature, respec tively, from the val ues reported 
in the discovery paper (|Skillen et al.l 12009?) is due pri- 
marily to the high-quality, high-resolution HIRES spec- 
trum analyzed. Moreover, the comparison of the results 
from the different stellar characterization methods, al- 
low us to derive uncertainties based on the internal er- 
rors, and including a systematic contribution, based on 
the range of derived stellar properties. Comparing the 
stellar spectroscopic properties derived in this paper to 
those in the literature, we find a hotter, and slightly 
more metal-rich star (T Bff = 5989 ± 51 K; [Fe/H] = 
0.06 ± 0.03 dex) than lSkillenetal] (f2009h . which is en- 
tir ely consistent with the stellar spectroscopic properties 
of iTorres et~aTI (|2012[) for WASP-13. The stellar mass 
that we derive (1.187 ± 0.065 M©) falls between the re- 
cent estimates for WASP-13 (1.09 ± 0.04 and 1.22 ± 0.1 2 
M©, respectively: iBarros et aHl2012t [Southworthll2012[ ). 
and above of that initially derived in the discovery pa- 
per (1.03 ± 0.10 M©). Correspondingly, the Saturn-mass 
planet WASP-13b is found to be also slightly larger and 
more massive (R p i = 1.407 ± 0.052 R Jup ; M pl = 0.500 ± 



0.037 M.jup ) with respect to the results from lSkillen et al.l 
(|2009j) . and lBarros et al.l (120121); but slig htly less massive 
and smaller than in iSouthworth] (|2012j ). Furthermore, 
we derive a younger planetary system (~4-5.5 Gyr) than 
before. While WASP-13 has certainly evolved off the 
zero-age main sequence, and might have not reached the 
overall contraction phase, it is equally likely that it has 
undergone contraction and has exhausted hydrogen in its 
core. This uncertainty in the stellar evolutionary stage 
of WASP-13 translates into an uncertainty in the stellar 
mass (see §3.3| and Fig. [3]) , which has not been typically 
accounted for. It is clear from our analysis that both a 
detailed stellar characterization, and a transit model in- 
cluding light curves and radial velocities are necessary to 
accurately determine the physical properties of planetary 
systems. 

It must be noted that there is very good agreement 
among the spectroscopic properties (Toff, [Fc/H], and 
log g) of WASP-13 derived from the four independent 
methods of stellar characterization for the unconstrained 
analysis. These stellar properties are also consistent 
with the independent measurement of log g constrained 
from the transit light curve (via p*), and of T e ff de- 
rived from the H a spectrum. Thus, this suggests that 
there is no significant systematic differences between the 
methods, and that unconstrained spectroscopic analyses 
give reliable stellar parameters for a star such as WASP- 
13. Although seemingly at odds with the conclusions 
from ITorres et al.l ()2012f ) , our methodology differs from 
their spectral synthesis-based analyses (see also £I3.2.1|) . 
Among these differences are the treatment of the micro- 
turbulence, the linelist, the sampling of a large parameter 
space in initial parameters, and the convergence crite- 
ria. Specifically our linelist includes the Na I D region, 
as well as the region of the Mg b triplet . It c ould be 
that the biases identified by ITorres et al.1 (|2012| ) due to 
the spectrocopic-log g on the other stellar parameters are 
not present in our implementation of SME. In addition, 
we have fixed the T G fj to the value derived from the H a 
analysis ( §3.1|) in the four stellar characterization meth- 
ods to assess the effect on the other spectroscopically- 
dctcrmincd stellar properties. In the case of WASP-13, 
we do not find any significant differences and the so- 
lutions are consistent with the unconstrained analyses. 
Although this is unsurprising given that all the derived 
temperatures are in agreement, the effect of fixing T c ff for 
different kinds of stars remains to be fully tested. With 
the larger HoSTS sample, a more robust conclusion, as 
to whether or not to constrain the spectroscopic anal- 
ysis of the planet hosts — with the stellar density from 
the transit light curve or with the T B s from a different 
temperature diagnostic — will be possible as a function 
of analysis method, quality of the data, and/or different 
stellar and planetary properties. 

Empirical trends between the physical proper- 
ties of planet hosting stars and their orbiting 
exoplanets have been previou sly identified and 

have been well stud ied (e.g. 

Fortnev fc Nettelmannl 120101: 



Baraffe et al.l 120101: iLaughlin et al.l 120111: 
2012t iFaedi et al.l 120111: iDemorv fc Seagerl 



iGuillot et alJ |200l 
iBouchv et aiT I201C 



Enoch et al.1 
20111) . 



For 

example, the relationships between planetary radius and 
stellar metallicity, as well as that between the planetary 
radius and the stellar irradiation, have been proposed to 
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Fig. 3. — Parameters of WASP-13 compared to the Yonsei-Yale stellar evolution models (Dcmarque ct al. 2001). The effective temperature 
and mean stellar density of WASP-13 is plotted as a solid circle. The six mass tracks have been interpolated in metallicity to [Fe/H]=+0.06 
to match the measured value of WASP-13. The mass of each track is labeled at the bottom of each line. Different evolutionary phases - 
main sequence evolution (solid black), overall contraction phase (dashed red) and post main sequenc e (blue do tted) — are given different 
line styles. Previou s measurement s of WASP-13 are a lso plotted in grey - diamond l|Barros et alj|2012f >. triangle (Southworth 2012), square 
IITorres et al.ll20l3) . and asterisk ISkillen et al.l[2009l ). The best mass for WASP-13 derived from these tracks is 1.245 Mq if the star is on 
the main sequence, but a 1.175 Mq star is equally likely if the star is slightly older and has already exhausted hydrogen in its core. 



probe planetary formation and st ructure. According t o 
theoretical models, like those of iFortnev et al.l (|2007|) . 
the planetary radius depends on the mass in the core of 
the planet. For example, a core-less Jupiter-mass planet 
that is dominated by its envelope has a radius that is 
several percent larger than a Jupiter-mass planet with 
a core- mass of a few tens of M®. At lower planetary 
masses, like in the Saturn-mass range (0.1 < M p \ < 0.5; 
lEnoch et al.ll2012|) . as the planets become dominated by 
the core-mass, the difference in radii for planets with 
and without cores seems more pronounced (see lower 
panel Fig. @| The higher stellar metallicity could lead to 
the formation of planet cores with more metal content, 
and t hus to smaller planetary radii (e.g., IGuillot et al.l 
I2006D . The planetary radius is also affected by the 
amount of irradiation received from the host star 
(iDemorv fc SeagerlrMTll lEnoch et al]|2012t iPerna et al.l 
120121) . Because of these different contributions, the 
effect of the stellar metallicity and irradiation on the 
formation and evolution of planets remains unclear. The 
structure of the planet and its environment need to be 
better constrained to be able to account for the diversity 



of physical properties of the known transiting systems. 
Thus, understanding the empirical relationships that 
have been previously identified between the stellar 
metallicity and stellar irradiation and the planetary 
radius may expand our knowledge on planetary systems. 

The top panel of Fig. @] shows the brightest (V < 14 
mag) transiting planets orbiting closest to the host star 
(Porb < 15 d) with planetary masses between 0.1 and 
12.5 Mj up , which compose the bulk of the HoSTS sam- 
ple. Doing a linear regression including the data uncer- 
tainties, we get a correlation coefficient of -0.24 between 
the stellar metallicity and the planetary radius. This 
known anticorrelation is more significant (with a corre- 
lation coefficient of -0.37) for the Saturn-mass planets 
shown in the lower panel as mentioned above. However, 
there is a strong correlation betw een the amount of s tel- 
lar flux received expressed by T eq ([Fressin et al.ll2007l see 
color-scale in lower panel) which has not been taken into 
account in the linear regression. In this paper, we do 
not attempt to characterize the observed trends except 
qualitatively. We note the evident improvement in the 
precision of the derived physical properties of the WASP- 



10 



Gomez Maqueo Chew et al. 



TABLE 3 

System parameters and la error limits derived from the MCMC 

ANALYSIS. 



Parameter 


Symbol 


Value 


Units 


Transit epoch 


T 


5305.62823 ± 0.00025 


dayst 


Orbital period 


P 


4.3530135 ± 0.0000027 


days 


Planet /star area ratio 


(R P i/R*) 2 


0.00844 ± 0.00016 




Transit duration 


t T 


0.1668 ±0.0010 


days 


Impact parameter 


b 


0.603 ± 0.025 




Stellar reflex velocity 


Ki 


0.0555 ± 0.0036 


km s — 1 


Center-of-mass velocity 


7 


9.8348 ± 0.0009 


km s _1 


Orbital eccentricity 


e 


0. 


fixed 


Orbital inclination 


i 


85.43 ±0.29 


degrees 


Stellar density 


P* 


0.306 ± 0.020 


P0 


Stellar mass 


M* 


1.187 ±0.065 


M Q 


Stellar radius 


h\ 


1.574 ±0.048 


R 


Orbital semi-major axis 


a 


0.0552 ±0.0010 


AU 


Planet radius 


R P i 


1.407 ±0.052 


R,iup 


Planet mass 


M pi 


0.500 ± 0.037 


M Jup 


Planet surface gravity 




2.764 ±0.038 


[cgs] 


Planet density 


Ppi 


0.180 ±0.020 


Pi 


Planet temperature 




1548 ± 22 


K 



t Given in BJD TDB - 2 450 000 as defined bv lEastman et al.l <201Cft . 



13 planetary system from our analysis (from gray-shaded 
area to fuchsia-shaded area). Our analysis of WASP- 13 
showcases how the HoSTS project will allow us to tighten 
the constraints rendered by the known transiting plan- 
ets, including the identification of any systematics, and 
correlations between the planet and stellar properties, 
as well as to assess the validity of the known empirical 
trends. 

We are continuing to acquire and analyze the high- 
quality echelle spectra and the long-slit spectra for T e g 



determination for the known transiting hosts. The data 
products of the HoSTS project are the derivation of a 
homogeneous set of stellar and planetary properties that 
will allow us to identify any biases in the parameters 
arising from the analysis methods and the quality of the 
data. This will enable us to significantly compare the 
physical properties allowing us to discover, derive, iden- 
tify trends among the planetary parameters exploring in 
more detail the planetary mass-radius relationship. 
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Fig. 4. — Top panel: R p i versus [Fe/H] for the bulk of the HoSTS sample taken from exoplanets.org (2012 Nov 19), and complemented 
with the literature. A linear regression to th e data gives a correlation coefficient of -0.24, which supports the previously observed trend 
fe.g.. iLaughlin et al. 2011; Enoch et al. 2012]). WASP-13 is marked by the filled circles with the shaded uncertainty areas: the grey box 
is for the [Fe/H] from Skillcn et al. (2009), and the from FSarros ct al] 1 120121) . and the fuchsia oval represents the properties derived 
in this paper. The uncertainties of the other data points are given by the grey crosses. Lower panel: We show the trend between R p \, 
stellar [Fe/H], and the planetary T cq (dependent on the stellar irradiation) for the HoSTS targets with planets in the Saturn-mass range 
(0.1 < M p i < 0.5 Mj up ). Although the two WASP-13 points overlap, the errors on the mctallicity are significantly smaller, showing the 
potential of the HoSTS project, by tightening the constraints rendered by the known transiting planets, including the identification of any 
systematics, and assess the validity of the observed trends. 
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